Search for CP- violation in Positronium Decay 
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CP-violation in the quark sector has been well established over the last decade, but has not been 
observed in the lepton sector. We search for CP-violating decay processes in positronium, using the 
angular correlation of {S ■ ki){S ■ ki x ^2), where S is the the positronium spin and fci, are the 
directions of the positronium decay photons. To a sensitivity of 2.2 x 10~'^, no CP- violation has 
been found, which is at the level of the CP-violation amplitude in the K meson. A 90% confidence 
interval of the CP-violation parameter (Ccp) was determined to be —0.0023 < Ccp < 0.0049. This 
result is a factor 7 more strict than that of the previous experiment. 



The material state of the universe results only be- 
cause of its matter-antimatter asymmetry. To under- 
stand the origin of this asymmetry is one of particle 
physics' most urgent goals. One cause could be the CP- 
violation (charge conjugation and parity) in the quark 
sector, which was found both in the K [l| and B [ij 
mesons. The CP- violating processes in B mesons have 
been measured precisely [3|] using its various decay pro- 
cesses at KEKB and PEP-II and they can be well de- 
scribed with the Kobayashi-Maskawa theory While 
understanding the origin of CP-violation in the quark 
sector is a great achievement, the resulting asymmetry 
may be too small to result in the amount of matter in 
the Universe All CP- violating observables in the 
Standard Model are proportional to the Jarlskog invari- 
ance, J= S2S3C1C2C3 sin(5) ^ 3 x 10~^, where Ci are Si 
are the mixing angles in the KM matrix, and S is the 
CP-violating phase. The CP-violating observables are 
also suppressed by the Higgs vacuum expectation value, 
since non-equilibrium can occur in the electro- 
weak phase transition in the Standard Model. It seems 
that these violations alone are too small to accout for 
the observed baryons in the universe. A new source of 
CP- violation or a new mechanism is necessary, and CP- 
violation in the lepton sector Q is a good candidate. 
Leptonic CP-violation has not yet been observed despite 
many studies at various neutrino factories 

Positronium (Ps) is a neutral system comprising an 
electron and a positron and is the analogue of the neu- 
tral meson. If CP is violated in the lepton sector, such 
neutral systems {e.g. Ps, muonium) are good test can- 
didates since the admixture of opposite CP eigenstates 
will occur. The search for CP-violation in Ps has been 
performed previously and no CP violation was found at 
the 1.5% level of uncertainty 0. 

If CP-symmetry is not conserved in Ps, the decay rate 
of the triplet Ps (o-Ps) is dependent on the CP- violating 
term. The number of o-Ps decay events would be ex- 
pressed as 

N ^ No[l + CcpiS-h)iS-ki X fc;)]exp(-t/T), (1) 



where S is the spin of the o-Ps, ki is the direction of 
the i-th energetic 7-ray from the decay (|A:i|=l), and 
T is the lifetime of o-Ps. Ccp is the amplitude of the 
CP-violation. The Standard Model prediction for Ccp 
is of the order of 10~^° [13] ■ This is due to photon- 
photon interactions in the final state. Another exotic 
decay contribution ha s been excluded and should be less 
than 2 X 10" 
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The discrete symmetry properties in the second term 
in Eq. ([T]) are as follows: All five vectors change sign un- 
der the T operation, while ki and S are a vector and an 
axial vector, respectively, under the P operation. Thus 
the angular correlation is CP-odd and T odd but con- 
serves CPT symmetry. This correlation would appear in 
general [l^l when the CP-violating effect exists in Ps. 

The analyzing power Q is defined as 



{S-ki){S-ki X fcs) 
P2 • sin 29 sin ip cos cj) 



(2) 



6 is the angle between the normal to the o-Ps decay plane 
and the spin quantization axis {S) of the o-Ps, is the an- 
gle between ki and k2 , and is the angle between fci and 
the projection of the spin quantization axis onto the o-Ps 
decay plane. Two parameters, 9 and -0, have been fixed 
during this measurement. CcpQifj)) "^iU be observed as 
an asymmetric function of if the CP-violation exists. 
P2 is the spin alignment of o-Ps (tensor polarization) and 
it can be expressed as 



P2 = 



iV+i +N0 + N.1 



(3) 



where Ni is the population of the iris ^ i Ps state. Usu- 
ally the spin alignment P2 is equal to zero since these 
three states are degenerate and thus their populations 
are expected to be the same, i.e. N^i = Nq = iV„i. A 
magnetic field is used to separate the triplet as follows: In 
a magnetic field, the o-Ps state with nis = and the p-Ps 
state mix, and the resultant states have different energies 
(Zeeman effect). This mixture depends on the strength 
of the magnetic field (0.2% for 5 kG). The lifetime of the 
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FIG. 1: Schematic diagrams of the experimental setup. The 
magnified view around the ^^Na source is shown in the circle. 
Left and right show the top and side view, respectively. 



resultant o-Ps state with mg = is greatly reduced to 
about 22 ns although the lifetime of the perturbed sin- 
glet is almost the same as that of p-Ps (125 ps) [l^l. On 
the other hand, the o-Ps states with nig = ±1 are not af- 
fected at all in a static magnetic field and their lifetimes 
arc still 142.05 ns in vacuum [l7, 18|. Since the = ±1 
and states contribute opposite alignment as shown in 
Eq. (jSj, we separate these Ps states on the basis of their 
different lifetimes in the external magnetic field. 

Figure [T] shows a schematic view of the experimental 
setup. A 1 MBq ^^Na positron source is sandwiched be- 
tween two sheets of thin plastic scintillator (NE102, t ~ 
0.1 mm). The emitted positrons pass through the scintil- 
lator and produce light pulses that are directed to photo- 
multipliers (2-inch fine-mesh HAMAMATSU H6614-70) 
by the light guide shown horizontally in both top and side 
views of Fig. [TJ The positrons form Ps when stopped in 
the targets of silica aerogel, 10 mm (diam.) x 5 mm 
(thickness), and density 0.11 g/cm'^. In order to reduce 
the pick-off annihilation of the Ps we use aerogel in which 
the surface of the primary grain is made hydrophobic, 
and wc also flush the aerogel with nitrogen gas (1 atm). 

The gamma-rays emitted from Ps decay are observed 
in four LYSO (Lui.gYo. 28105) crystals (diameter = 30 
mm and length = 30 mm) shown as the hatched boxes 
in the top view. The four detectors are placed such as 
to make three pairs with opening angles equal to 150°. 
These three pairs are used to detect the 37 decay of o-Ps. 
The scintillation pulses of the LYSO crystals are detected 
with fi-inch fine- mesh photomultipliers (HAMAMATSU 
H8409-70). The average energy resolution of the LYSO 
detectors is 60 keV (FWHM) at 511 keV, and a good 
timing resolution of 1.2 ns (FWHM) is obtained for 511 
keV gamma-rays. Since the mean atomic number (Z) of 
the LYSO crystal is large, the photoelectric process dom- 
inates compared to the Compton scatter process. This is 
an advantage for tagging fci and k2 correctly. 



A magnetic field to perturb the Ps spin states is gener- 
ated using two 80 mm (diam.) x 50 mm (thickness) Nd 
permanent magnets with iron cones and yokes as shown 
in the side view. The magnetic field is 4.9 kG and is di- 
rected at 30° from the direction of the normal to the 37 
decay plane. Its uniformity is about 10% over the volume 
of the silica aerogel. The observed unperturbed and per- 
turbed o-Ps lifetimes are 126 ns and 22.5 ns, respectively. 
These lifetimes are consistent with the expected lifetimes 
in 1 atm nitrogen gas and a magnetic field of 4.9 kG. 

As shown in the top view, all scintillators are located 
on a turntable which turns with respect to the magnetic 
field. Using this turntable, we can vary </> in Eq. Q inde- 
pendently and measure the (/(-dependence of the number 
of events. The turntable is rotated by a stepping motor 
and its rotation angle is controlled using a CAMAC sys- 
tem. The angle of the turntable is measured by an optical 
sensor with 0.2° accuracy. The experimental apparatus 
is carefully assembled to avoid geometrical asymmetry 
and the assembly accuracy is about 0.5 mm. 

In order to select Ps decay events data acquisition logic 
is set up as follows: When at least two signals from the 
LYSO scintillators arc coincident within 20 ns, and then 
when this coincidence is within —100 ns to 600 ns of 
the timing of the plastic scintillators, data aqcuisition is 
triggered. The coincidence between the LYSO detectors 
reduces accidental counts due to the radioactive decay 
of in the crystal. A charge ADC (REPIC RPC- 

022) is used to measure the energy information of the 
plastic scintillators. The outputs of the LYSO detectors 
are recorded with another charge ADC (CAEN C1205). 
The time information of the plastic and LYSO scintilla- 
tors is recorded using a direct clock (2 GHz) count type 
TDC [13. 

Five runs with different alignments of the LYSO de- 
tectors were performed to check systematic effects de- 
pendent on the position of the LYSO detectors. The 
total data acquisition period was about six months. The 
trigger rates were about 1.3 kHz. During the data acqui- 
sition, the turntable was rotated around 30° every hour. 
Energy and time calibrations were also performed every 
hour. The 7-ray peak at 511 keV and the zero energy 
peak were used to calibrate the LYSO detectors. The 
room temperature was maintained within 22 ± 1°C in 
order to maintain good stability during the data acquisi- 
tion. 

Figure [5] shows the 2-dimensional distribution of 7- en- 
ergies obtained with the 150°-pairs of the LYSO detec- 
tors. 7-rays from 37 decays are distributed in the arc 
whose shape is determined by the opening angle of the 
detector pairs. To select 37 events, the following energy- 
window is applied: 

680 [keV] <Ei+E2< 920 [keV], (4) 

280 [keV] <E2<Ei<5n [keV], (5) 

where Ei is energy of i-th energetic 7-ray emitted from 
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FIG. 2: 2D energy distribution of o-Ps decay 7's obtained by 
150° pairs of the LYSO detectors. 



FIG. 4: Measured asymmetry averaged over the detector pairs 
and five runs. The fine shows the expected angular depen- 
dence if the CP violating coefficient Ccp = 0.01. 
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FIG. 3: Time difference between the plastic scintiUators and 
the coincidence signals of the LYSO scintillators after the en- 
ergy selection is applied. The figure in the right corner is the 
magnified view of the sharp peak of the prompt annihilation. 



the o-Ps decay. This selection is superimposed in Fig. [2] 
Figure [3] shows the time difference between the plas- 
tic scintillator signal and the coincidence signals of the 
LYSO detectors after the energy selection is applied. A 
sharp peak from prompt annihilation is followed by the 
exponential curves of o-Ps and then the constant acciden- 
tal spectrum. The figure in the right corner shows the 
magnified view of the sharp peak of the prompt annihila- 
tion. A good timing resolution of ct = 0.9 ns is obtained. 
A time window from 50 ns to 130 ns is required to en- 
hance the ma = ±1 Ps states and a spin alignment P2 ~ 
0.87 is obtained. 

The asymmetry, is calculated using the numbers 

of events, N((j)), from which the accidental contribution 
estimated with the time spectrum shown in Fig. |3] is 



N{(j>) ~ N{(l)+ 180°) 
N{(t)) + N{(t) + 180°)' 



(6) 



This asymmetry is equal to CcpQ{(f>) since the numbers 
of events are expressed as Eq. ^ and the analyzing 
power Q{(j)) depends on coscj) as shown in Eq. ([2]). The 
differences between the LYSO detector pairs do not af- 
fect the measured asymmetry since an independent asym- 
metry measurement is performed for each detector pair 
(thanks to the rotating table). The effect of Compton 
scatter by the magnet near the ^^Na source is obviously 
the same between N{(j)) and N{(l)+ 180°) since the mag- 
net system is symmetric with respect to the origin of the 
experimental setup. In addition, averaging the asym- 
metries of the same analyzing power over the pairs of 
the LYSO detectors cancels out systematic effects due to 
other unknown sources. In Fig. |4l the asymmetries are 
plotted as a function of (p. If CP-symmetry is violated, 
the A{(j)) distribution would be a cosine curve (superim- 
posed in the figure) whose amplitude is proportional to 
the magnitude of the CP- violation. 

The obtained data is consistent with the flat distribu- 
tion as shown in Fig. |4] and thus no CP- violating effect 
is observed. The Ccp coefficient is obtained by dividing 
the measured asymmetries by the analyzing power. The 
analyzing power shown in Eq. (2) is the ideal case, but 
the real analyzing power of this setup is affected by var- 
ious factors: the finite solid angle (—15% with respect 
to the ideal case), the energy resolution (—2%) of the 
LYSO crystals, and the finite size of the Ps creation po- 
sition (—1%). These were estimated using a Monte Carlo 
simulation. The angular correlation of the 7-rays from 
the nis = 0, ±1 states of o-Ps is also taken into account 
using a theoretical calculation [l^. The obtained coef- 
ficient averaged over (j) is Ccp = 0.0013 ± 0.0021(stat.), 
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TABLE I: Summary of the systematic errors. 



Source 


Systematic error 


Table angle accuracy 


±0.00039 


Center alignment 


±0.00025 


Non-uniformity of aerogel / magnetic field 


±0.00011 


Decrease of Z?"*" decay rate 


±0.00030 



and Ccp is consistent with zero within the statistical 
uncertainty. 

Estimates of various systematic errors are summarized 
in Table HI The largest contribution is the uncertainty 
in the angle of the turntable. The turntable angle is 
directly related to 0, and the effect of Compton scat- 
ter by the magnet can not be canceled out if the angle 
between N{(j)) and A^(0 + 180°) differs from 180°. The 
number of events differed by ±10% when the table was 
rotated by 90° due to this geometrical effect and this 
difference is applied to the 180°-case conservatively. An 
uncertainty (0.2°) in the table angle could contribute an 
uncertainty of about 2 x 10~^ to the event number. This 
error propageted to Ccp is 3.9 x 10~^. Another predom- 
inant systematic error is produced if the center of the 
turntable and the magnetic field are misaligned. In such 
cases, the number of events which pass the time win- 
dow changes since the mean magnetic field over the Ps 
creation area depends on the table angle. From the mea- 
sured magnetic field, and the assembly accuracy of the 
experimental setup (0.5 mm), we can estimate the uncer- 
tainty of the mean magnetic field, and it contributes to 
the error on Ccp by 2.5 x 10~^. Non-uniformities in the 
scilica aerogel and the magnetic field would also produce 
asymmetries, since they also make the effective strength 
of the magnetic field depend on the table angle. The 
last predominant source is the decrease in the decay rate 
of ^^Na {Ti/2 = 2.6 yr). In this experiment, N[4>) and 
N{<f) + 180°) are recorded at different times (time differ- 
ence = 6 h). The decay rate of ^^Na decreases slightly 
over this time (2.6 x 10"^). This results in a small de- 
crease in the number of events and causes an asymmetry. 

The systematic errors discussed above are regarded as 
independent contributions such that the total systematic 
error can be calculated as their quadratic sum. The com- 
bined result with the systematic errors is 

Ccp ^ 0.0013 ± 0.0021(stat.) ± 0.0006(sys.). (7) 

Cc p is consistent with zero and a 90% confidence interval 
of 

- 0.0023 < Ccp < 0.0049, (8) 

is obtained. This result is 7 times more strict than that of 
the previous experiment Q and is of a similar sensitivity 
to that of the Kaon CP-violation. 



In order to perform more accurate experiment you need 
more statistics but you can not simply use a radioactive 
source with higher intensities because it causes pulse pile- 
up in detector systems. A possible way to increase statis- 
tics is to use detector arrays with nearly 47r sr of coverage. 
In addition, as for systematics, an electromagnet can be 
useful because it can also cancel many systematic effects 
by reversing its polarity. 

A CP-violating decay process in positronium has been 
searched for using the angular correlation of (S ■ ki){S ■ 
ki X k2), where S is the the positronium spin and fci, k2 
are the directions of the decay photons. LYSO scintilla- 
tors are used to measure the photon energies and these 
are set up to turn with respect to the magnetic field in 
order to cancel out various systematic errors. To a sensi- 
tivity of 2.2 X 10~^, no CP-violation has been observed. 
This level is similar to the CP- violation amplitude of the 
K meson. 
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